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1.0  AIM  System  Concept 


Strategic  defense  of  the  continental  United  States  (CONUS)  depends  upon  our  ability  to  detect  and 
track  potential  threats  approaching  at  long  ranges,  allowing  adequate  time  for  assessment  and  reac¬ 
tion.  To  address  this  critical  need  for  early  warning  surveillance,  well  beyond  the  line  of  sight 
(LOS)  of  traditional  radars,  the  U.S.  Air  Force  has  developed  and  is  currently  operating  the  Over- 
the-Horizon  Backscatter  (OTH-B)  radar  surveillance  system  (AN/FPS-i  18)  [1],  In  addition  to 
OTH-B,  the  U.S.  Navy  has  developed  the  Relocatable  Over-the-Horizo,i  Radar  (ROTHR),  the 
AN/TPS-7 1  [2,3].  Operational  results  from  the  initial  east  coast  facility  indicate  that  OTH-B  will 
serve  as  a  valuable  cornerstone  to  the  early  warning  surveillance  capability  required  for  CONUS 
defense.  The  threats  of  particular  cunent  interest  are  penetrators  that  exhibit  small  radar  cross 
sections,  including  air  launched  and  submarine  launched  cruise  missiles  (ALCM  and  SLCM).  This 
threat  will  become  more  important  in  theater  warfare  as  well  over  the  next  few  decades.  The  ability 
to  supplement  current  wide  area  surveillance  with  reliable,  all  weather,  continuous  detection  and 
tracking  of  low  flying,  low  observable  targets  is  paramount  to  the  successful  implementation  of 
both  strategic  and  tactical  air  defense  strategies. 

Recognition  of  the  performance  limitations  associated  with  the  natural  ionosphere  motivated  a 
number  of  investigators  (eg.  A.  Drobot  and  D.  Papadopoulos;  A.  V.  Gurevich  [4,5])  to  propose 
the  creation  of  an  Artificial  Ionospheric  Mirror  (AIM)  in  the  upper  atmosphere,  in  order  to  reflect 
ground-based  radar  signals  for  OTH  surveillance.  The  AIM  is  produced  by  beaming  sufficient 
electromagnetic  power  to  the  lower  ionosphere  (around  70  km)  to  enhance  the  in  situ  ionization 
level  to  107-  10K  electrons/cm3,  thereby  providing  an  ioniced  layer  capable  of  reflecting  radar 
frequencies  of  5  -  90  MHz.  This  paper  presents  a  baseline  AIM  system  concept  and  an  associated 
performance  evaluation,  based  upon  the  relevant  ionization  and  propagation  physics  in  the  context 
of  air  surveillance  for  the  cruise  missile  threat.  Results  of  the  subject  study  indicate  that  a  system 
iiNing  this  concept  would  both  complement  and  enhance  the  performance  of  the  existing  skywave 
OTH  radars,  as  illustrated  in  Figure  1.  Performance  analyse:,  for  the  projected  system  show  that  it 
has  the  pc  tential  for  providing  reliable  and  consistent  detection  and  tracking  of  the  cruise  missile 
threat. 

Tire  performance  characteristics  of  a  stand-alone  AIM  system  are  highlighted  by 

•  200- 1 200  km  detect  ion/tracking  ranges  (depending  on  target  altitude), 

•  azimuth  coverage  over  360°,  or  any  fractional  sector, 

•  beam  grazing  angles  of  less  than  10  degrees  (for  65-80  km  high  AIM), 

*  good  horizontal  polarization  control, 

*  operation  at  RFs  in  the  HF  to  low  VHP'  (to  about  90  MHz),  and 

*  90%  detection  probability  of  a  -25  dBsm  target  at  1 .000  km  range,  providing  30  dli 
of  margin  for  typical  iow  observable  threat  projections  at  VHP. 


The  performance  of  a  normal  skywave  OTH  system  can  be  enhanced  by  providing  an  additional 
AIM  sector  (Figure  2),  which  has  sustained  operation  in  the  20-30  MHz,  regime  independent  of 
time  of  day.  latitude,  look  angle,  and  ionospheric  state.  Tire  introduction  of  an  AIM,  located  at  the 
appropriate  altitude,  will  create  a  reliable  and  predictable  reflection  of  the  HF  energy  between  the 
radar  site  and  the  area  of  it  tercst.  An  AIM  adjunct  to  a  conventional  HF  skywave  radar  can  offer 
substantially  improved  performance  with  regard  to: 

•  filling  in  the  range  hole  that  exists  out  to  about  1000  km  due  to  the  minimum  1  IP 
hop  distance. 

♦  mitigation  of  auroral  effects  in  polar  directed  surveillance  sectors. 


I 


minimum  grazing  angle 


WITHOUT  AIM  WITH  AIM 


Blinded  by  ionospheric  vagaries 

Skip  zone  permits  undetected  nearby  launch 

Small  targets  unseen  at  night 

Large  clutter  back  scatter 


Independent  of  ionospheric  vagaries 
Skip  zone  filled  in 
Nightime  performance  maintained 
Horizontal  polarization  reduces  clutter 


Figure  1.  An  AIM-Based  Radar  Will  Serve  as  a  Valuable  Complement  to  a  Skywave  OTH  Radar 

•  sustained  operation  through  periods  of  increased  sunspot  activity  and  other  iono¬ 
spheric  degradations, 

•  availability  of  the  upper  end  of  the  HF  spectrum  during  the  diurnal  ionospheric 
cycle,  and 

•  improved  detection  of  LO  targets  through  frequency  selection  and  positive  polar¬ 
ization  control. 

In  addition,  an  AIM  based  system  is  not  restricted  to  the  HF  band,  but  can  operate  in  the  lower 
VHF  band,  which  lias  several  advantages  over  HF: 

•  less  crowded  band  permits  broader  bandwidths,  resulting  in  better  resolution  and 
increased  sensitivity. 

•  lower  noise  temperature,  allows  higher  signal-to-noise  ratios  and  improved 
detection  performance. 

•  resonant  frequencies  for  cruise  nessilc  types  of  targets,  allowing  long  range 
detection  of  otherwise  very  iiT'.mlt  to  detect  threats,  as  illustrated  in  Figure  2. 

Specific  issues  influencing  an  AIM  system  design  fall  into  one  of  three  primary  categories: 

•  AIM  heater  design  trade-offs,  which  determine  how  well  one  can  control  the 
production  of  an  AIM  layer, 

•  AIM  RF  propagation  effects,  which  allow  one  to  assess  the  quality  of  the  reflected 
radar  wave  (eg.  absorption,  wavefront  distortions,  dopplcr  spreading,  and  Faraday 
rotation),  and 

•  Environmental  and  operational  requirements,  which  specify  the  threat  and  mission 
context  in  which  the  AIM  system  must  perform. 

Ibis  report  will  present  results  concerning  each  of  these  issues,  highlighting  and  quantifying  those 
mosf  critical  in  determining  the  bottom-line  radar  system  performance. 
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Figure  2.  AIM  Provides  an  Additional  High  Resolution,  Highly  Available  Surveillance  Sector 


3 


2.0  AIM  Formation  and  Control 


2.1  Background 

Results  from  theoretical  analysis  and  numerical  modelling  indicate  that  a  useful  AIM  consists  of  an 
ionized  layer,  with  an  electron  density  gradient  increasing  with  altitude  from  the  ambient  state  to  a 
density  with  plasma  frequency  corresponding  to  the  radar  frequency.  The  AIM  layer  is  oriented 
with  an  inclination  (~  40°-45°)  and  has  a  slight  curvature  to  allow  the  radar  wave  to  refract  over  a 
15°  elevation  beamwidth,  providing  range  coverage  from  200  to  1,200  km.  In  addition,  the  AIM 
azimuth  orientation  is  incrementally  stepped  over  the  desired  angle  of  surveillance. 

Figure  3  illustrates  the  system  timing  dependence  between  the  AIM  heater  operation  and  the 
surveillance  sequence  of  the  radar.  The  illustration  shows  the  time  coordination  of  the  AIM  heater 
and  the  radar  at  three  time  scales.  The  coarsest  scale  (bottom)  indicates  the  sequence  of  patch 
creation,  followed  by  radar  dwell.  Using  this  strategy,  the  direction  of  the  radar  dwell  can  be 
changed  at  each  dwell  (by  creating  a  new  patch  orientation),  thereby  achieving  a  scanning  radar 
beam.  The  middle  time-line  shows  a  single  AIM  creation/dwell  period,  in  which  the  AIM  is 
created  (approximately  10  to  50  ms),  followed  by  a  short  cooling  time  (less  than  .1  ms),  and  a 
radar  dwell  for  the  life  of  the  cloud  (~  Is).  The  size  and  orientation  of  the  AIM  relative  to  the  spot 
size  of  the  focussed  RF  heater  beam  (-36  tn  diameter  for  this  case)  requires  that  the  AIM  creation 
tie  performed  with  a  raster  scanning  procedure.  Figure  3  illustrates  the  time  sequence  for  this 
scanning  procedure,  in  which  a  complete  AIM  is  formed  by  moving  the  heater  focus  along  a 
stepped  sequence  of  horizontal  sweeps  (each  sweep  produces  a  bar  of  ionization  the  width  of  the 

patch  and  takes  -  70  to  350  ps).  The  altitude  at  which  ionization  occurs  is  determined  by  where 
the  power  density  exceeds  the  breakdown  threshold.  Positive  control  of  this  altitude  is  achieved  by 
using  a  large,  partially  filled,  array  and  focussing  the  beam  so  that  the  AIM  is  formed  on  the  front 
side  of  the  focus.  The  focal  position  is  controlled  so  as  to  follow  a  surface  contour  having  the 
desired  azimuth  and  elevation  orientation  and  a  slight  curvature  to  allow  for  radar  beam  divergence 
over  the  azimuth  elevation  required  for  surveillance,  litis  procedure  is  referred  to  as  "painting" 
ami  is  nolionally  illustrated  in  Figure  4.  Based  on  our  current  understanding  of  ionization  physics 
and  t\a .  „»t  technology,  a  typical  AIM  will  require  on  die  order  of  a  10  to  50  ms  to  paint. 

A  central  issue  concerned  with  evaluating  the  feasibility  of  the  AIM  concept  is  how  reliably  one  can 
create  the  AIM  within  tolerances  necessary  for  useful  reflection  of  the  radar  signal.  Issues  that 
directly  impact  usefulness  include:  size,  shape,  orientation,  uniformity,  smoothness,  peak  electron 
density,  steepness  of  the  density  gradient,  and  density  lifetime.  In  order  to  adequately  address 
these  issues,  one  needs  a  comprehensive  understanding  of  the  phenomenology  (validated 
experimentally)  and  a  careful  analysis  of  the  lieater  system  design  parameters  as  they  impact  plasma 
formation  control.  The  underlying  physics  has  been  a  major  focus  of  study  and  the  results  are 
reported  in  (6),  The  system  trade-off  considerations  indicate  that  the  -MM  can  be  created  using 
current  technology,  with  the  necessary  radar  reflective  characteristics  to  provide  wide  area 
surveillance  and  early  detection  and  tracking  of  the  A1.CM/S1.CM  threat  of  the  future.  Specific 
results  of  the  system  trade-off  studies  are  presented  in  the  following  discussions. 

2.2  Vertical  Electron  Density  Profile  Control 

As  indicated  abos-c,  altitude  control  of  the  AIM  formation  is  achieved  by  utilizing  a  large  t.xrussed 
heater  array,  so  that  ionization  hegms  on  die  front  edge  of  the  focus  (power  density  is  increasing 
with  altitude).  Breakdown  initiates  where  the  power  density  reaches  a  threshold  level  (30  kW/nr 
#  69  km  altitude  and  heater  frequency,  fh  425  MHz),  and  rapidly  increases  in  ionization  rate 
above  this  level.  As  the  number  of  electrons  increase,  absotption  of  the  healer  wave  begins  to 
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Figure  3.  AIM  System  Timing  Diagram 
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Figure  4.  The  Large.  Sparse  Heater  Can  Paint  a  Snwoth.  Homogeneous  AIM  Patch 
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reduce  the  power  reaching  altitudes  above  the  critical  altitude  and  continues  reducing  the  power 
level  until  it  drops  below  the  threshold.  This  "self-absorption"  mechanism  limits  ionization  levels 
to  plasma  frequencies  well  below  the  heater  frequency.  Figure  5  illustrates  a  time  evolution  of  the 
ionization  process  as  a  function  of  altitude.  TTiis  figure  clearly  shows  the  limiting  of  the  peak 
density  and  the  steepening  of  the  electron  density  profile  over  time.  This  steepening  is  due  to  the 
fact  that  self-absorption  causes  the  ionization  to  stop  at  the  higher  altitudes  first  and  then  to  walk 
back  toward  the  critical  altitude.  Upon  examination  of  the  graph,  we  note  that  the  critical  power 
density  for  this  situation  is  approximately  45  dBW/m2.  Moreover,  when  the  power  density  has 
attenuated  below  42  dBW/m2  at  70  km  altitude,  there  is  no  further  ionization  at  that  altitude. 


While  many  factors  impact  the  resulting  electron  density  profile,  the  major  influencing  factors, 
within  control  of  the  system  designer,  are  the  heater  frequency,  dwell  time,  and  the  power  density 
gradient.  Increasing  the  heater  frequency  reduces  the  absorption  rate,  thereby  allowing  higher  peak 
plasma  densities.  Longer  dwell  time  allows  the  the  profile  to  walk  further  back  toward  the  critical 
altitude,  hence  increasing  the  final  density  gradient.  The  third  factor,  power  density  gradient, 
determines  how  rapidly  the  power  density  is  increasing  beyond  the  threshold  density  and 
consequently  how  much  absorption  is  required  to  reduce  the  field  strength  down  to  the  critical 
level.  Wlide  higher  power  density  gradients  produce  higher  peak  plasma  frequencies,  increased 
gradients  are  obtained  at  the  expense  of  smaller  beam  widths  and  thus  smaller  ionization  areas. 
Approximate  relationships  between  these  factors  and  the  electron  density  profile  parameters  have 
been  empirically  determined  for  AiM  formation  at  70  km  altitude  (71.  The  peak  electron  density  is 


Power  Density  (dBw/m  2) 

42  44  46  48  50  52 


(IJ 


Ne  =  min  j  Neol.5  exfa.4  x  10  i2(pf/fif  r),  2  x  10 1  fhA  VP, 15 


and  the  electron  density  gradient  is 

14  exp  |6.5  x  10  6  f58  7.5  x  104/*  35  max  (v/%,4]  j ,  (2) 
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power  density  at  the  power  peak  (W/m  ) 
heater  frequency  (GHz) 

pulse  length  (Msec) 

ambient  electron  density  (cm‘3)  =  102  cm 3  at  70  km 
power  density  gradient  at  the  electron  density  peak  (dB/km). 

In  each  case  the  first  term  in  min  { •,  •  }  represents  the  growth  phase  of  the  cloud,  while  the  second 

term  represents  the  maximum  or  clamped  state  of  the  cloud.  As  t  increases,  the  growth  formula 
approaches  and  finally  reaches  the  value  of  the  damped  formula,  at  which  point  the  second  term 
becomes  the  valid  expression.  Since  these  relationships  are  only  valid  for  a  neutral  density 
corresponding  to  70  km  altitude,  in  order  to  support  a  full  system  design  trade-off,  they  will  be 
expanded  to  account  for  altitude  dependence  (ie.  sensitivity  to  neutral  density). 

As  will  be  discussed  in  Section  3.0,  these  profile  characteristics  have  important  implications 
relative  to  utilizing  the  AIM  as  an  RF  reflector,  determining  both  the  maximum  useable  frequency 
(MUF)  and  the  amount  of  absorption  and  wavefront  distortion  incurred  while  the  wave  transits  the 
cloud  Therefore,  the  AIM  system  engineer  must  carefully  trade-off  the  relevant  heater 
characteristics  and  resulting  performance  degradations  against  overall  cost  considerations.  In  order 
to  adequately  perform  such  a  trade-off,  the  relation  between  heater/AIM  characteristics  and  system 
performance  must  be  quantified.  Quantification  of  two  major  factors  is  discussed  in  the  following 
two  subsections,  while  the  remaining  issues  are  addressed  in  Section  3.0. 

2.3  Healer  Control  Errors 

As  mentioned  above,  formation  of  the  AIM  plasma  will  be  achieved  by  RF  radiation  from  a  very- 
large  ground-based  phased  array.  Since  neutral  density  atm  electron  density  spatial  variations  ate 
less  than  05*£  over  AIM  scale  lengths  JSJ.  the  primary  source  of  plasma  irregularities  will  be  due 
to  heater  scanning  and  control.  The  objective  of  the  follow  ing  analysts  and  associated  discussion 
is  to  derive  and  evaluate  some  rough  measures  of  sensitivity  between  the  inlterent  uncertainties  in 
generating  the  individual  heater  element  waveforms,  the  resulting  perturbations  it-  AIM  breakdown 
location,  and  subsequent  impact  mi  radar  system  performance.  Deviations  in  waveform  generation 
at  the  element  level  result  in  corresponding  deviations  in  the  desired  power  density  profile  in  the 
breakdown  region.  Three  specific  potential  sources  of  error  are  considered: 

•  simple  additive  muse. 

•  amplitude  distortions,  and 

•  phase  distortions. 

While  trie  following  discussion  specifically  addresses  issues  related  to  3  hrge  digitally  controlled 
solid  state  heater,  the  results  can  be  readily  applied  to  high  poorer  microwave  sources. 

Examples  of  potential  causes  of  these  errors  include  quantization  errors  in  the  digital  circuitry, 
dsgital-to-analog  conversion  errors,  power  amplifier  gain  deviations  and  distortions,  and  timing 
emus  in  the  waveform  control  circuitry.  Each  of  these  factors  ir  considered  relative  to  what  is 
readily  achievable  w  ub  current  technology  and  how  these  performance  limitations  impact  the  patch 
integrity. 


dNe 

dz 


=  min 
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Results  of  this  analysis  indicate  that  although  the  power  density/breakdown  dependencies  require 
relatively  accurate  control  of  the  power  density  profile  (c.g.  to  within  .003  dB,  which  corresponds 
to  a  SNR  at  the  focal  point  of  68  dB),  the  large  number  of  independent  heater  elements  resuk  in 
significant  integration  gains  (e.g.  a  configuration  of  40,000  elements  reduces  the  68  dB 
requirement  at  the  focal  point  by  46  dB  to  only  22  dB  at  the  element  level).  A  summary  of  e 
impact  of  this  integration  gain  on  what  is  achievable  with  currently  available  hardware  technol  / 
is  summarized  in  Table  1. 

Table  1.  Heater  Error  Budget  Achievable  with  Currently  Available  Hardware 


Error  Factor 

Deviation  in  AIM  Location 

Quantization  Noise 

0.13  m 

D/A  Noise 

.009  m 

Power  Amplifier  Distortions 

<  ,13  ni 

Phase  Deviations 

,33  m 

36  m 

For  the  surveillance  radar  frequency  operating  at  50  Mhz,  the  patch  irregularities  indicated  by  the 
cumulative  error  is  on  the  order  of  .06  wavelengths.  As  indicated  by  Greene’s  analysis  [9],  an 
aperture  with  rms  phase  errors  of  this  order  of  magnitude  will  experience  minimal  performance 
degradation  (less  than  .3  dB  reduction  in  gain,  less  than  .2%  increase  in  beamwidth,  and  root- 
mean-square  (rms)  beam  canting  less  than  10%  of  the  beamwidth).  Furthermore,  when  path 
averaging  of  the  radar  wave  within  the  plasma  is  taken  into  account,  tlie  results  indicate  that  the  net 
effect  on  radar  performance  is  negligible  (~  .02  dB). 

In  the  following  discussion,  we  derive  the  relationship  between  clement  waveform  errors  and  ’.he 
resulting  error  in  plasma  breakdown  location.  Each  potential  error  source  is  related  to  an 
equivalent  dement  waveform  error.  These  two  relationships  are  then  combined  to  determine  the 
sensitivity  of  the  AIM  breakdown  location  to  each  noise  factor.  The  consequent  AIM  location 
emirs  are  used  in  Section  3.0  to  evaluate  the  impact  these  factors  have  on  radar  detection. 


2.3.1  Sensitivity  Relationships 

Figure  6  illustrates  tire  basic  patch- hearer  geometry  with  tire  indicated  parameters  defined  as 


<Sn 

£(U) 

£*(f) 

& 


desired  location  of  the  array  focus,  relative  to  the  center  of  the  heater, 
location  of  the  nth  heater  clement,  relative  to  the  center  of  the  heater, 
desired  electric  field  at  location  r  and  lime  i. 
desired  clement  wavefonn. 

dement  amplitude  gain  pattern,  and 

It  -  wwt  direction  vector  kom  tire  nth  element. 


For  a  heater  array  with  N  elements  and  c  denoting  the  speed  of  light,  tire  ideal  dcctnc  field  intensity 
is  given  by  ihe  summed  contributions  from  each  clement. 

fc-fco  .  ijoJ®* 


I  M 


*  £  s(sj£Jt  hr  j^cj  . 


K 


Figure  6.  Geometry  of  Heater  Array  Sensitivity  Analysis 

where  P  is  the  polarization  orientation  of  the  electric  field  and  Ifcll »  Hill. 

Equation  (3)  provides  the  relationship  between  the  electric  field  intensity  at  location  r  and  the  ideal 
element  waveforms.  For  a  non-ideal  waveform,  we  assume  that  g(j^)En(t)  has  been  perturbed  by 
an  additive  noise  process,  rjn,  which  has  zero  mean  and  variance  £(  irjn(r)l2)  =  cr^2,  where  £{  •  ) 
denotes  the  statistical  expectation.  The  actual  contribution  to  I&,t)  by  each  element  is 

g{<pn)Eti)+Ut). 

In  the  analysis  to  follow,  we  assume  the  element  errors  are  uncorrclated  from  element  to  element. 
This  appears  to  be  a  reasonable  assumption  for  the  error  sources  mentioned  above.1  Given  this 
assumption,  the  desired  power  density  and  potential  noise  power  at  the  focal  point  of  the  heater  arc 

Pf=t{\Eiif,t)\  2)  =  ^NzE2  (4) 

and 

N  2 

cj/= £  -Na\.  (5) 

U*l 

In  equations  (4)  and  (5)  we  have  approximated  the  1/Hc  -  inll2  spreading  loss  to  be  the  same  for 
each  clement,  thus  omitting  the  common  l/(4jt/?2)  factors  from  (4)  and  (5).  In  addition,  to 

1  One  possible  exception  is  the  timing  circuitry.  Timing  errors  due  to  system  dock  jitter  would  correlate  across 
all  of  the  heater  elements.  This  is  discussed  in  more  detail  at  a  later  point  in  the  paper. 
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simplify  the  analysis,  we  assume  that  En  is  generated  such  that  the  magnitude  of  the  g{^n)En  are 

equal  (\g(Jhn)En\  =  E).  Combining  (4)  and  (5)  gives  the  signal-to-noise  ratio  (SNR)  of  the  power 
density  at  the  focal  point  of  the  heater 

M£l2 

SNR  =  =  N  SNRe ,  (6) 

2  O? 

where  SNRe  is  the  signal-to-noise  ratio  at  the  output  of  each  individual  heater  element. 


In  order  to  relate  the  SNR  at  the  focus  of  the  array  to  resulting  irregularities  in  the  patch  breakdown 
location,  we  will  consider  the  effects  of  deviations  from  the  desired  power  density  profile  on 
breakdown  location.  Figure  7  illustrates  a  representative  power  density  profile  in  one  dimension 
with  the  corresponding  plasma  electron  density  after  breakdown.  In  the  region  of  breakdown,  the 
power  density  is  approximated  by 


\Qlog\dP)  =  az  +  J 3, 


(7) 


where  z  is  the  altitude  from  the  heater  array.  Using  (7),  a  small  deviation  in  power,  A P,  will  result 

i>'  i  ..’ritude  offset,  Az;  representing  the  difference  between  the  desired  location,  at  which  the 
pc  wev  aensity  reaches  the  critical  breakdown  level,  Pc,  anti  the  actual  location  of  the  critical  powei 
density, 

Az  =  &togl  o(l+^).  (8) 

U  f  c 


A  linear  approximation  in  (8)  gives 


Az 


_ _ AE. 

atn(  10)  Pc  ‘ 


(9) 


If  the  r/n(t)  's  are  independent  random  processes,  the  central  limit  theorem  implies  that  Az  is  a  non- 
central  chi-squared  random  variable,  with  mean  and  variance  derived  from  (9) 


eU->! 


and 


_ ULL,. 

aln(  10) 


.ol 

Pc 


oi 


A? 


1 4.  <£. 

2Pj 

a  /«(!()) 

Pc' 

a/fl(IO) 

SNRC' 


(10) 

(il) 


In  ( 10)  and  (II).  represents  die  variance  in  the  combined  field  error  at  the  breakdown  altitude, 
due  to  all  elements  and  SNRC  is  the  signal  to  noise  ratio  at  the  breakdown  point.  SNRC  will  be  less 
than  at  the  focus,  because  dcfocussing  of  the  elements  causes  the  expected  field  intensity  to  be 
reduced,  while  the  sum  of  independent  clement  errors  is  not  affected  by  the  phase  difference 
between  elements.  Combining  equations  (6),  (10),  and  (11)  gives  the  variance  in  the  error 

between  the  average  breakdown  location,  z( ,  and  tlic  actual  location.  zc  +  Az, 
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- 20 - ,  (12) 

aln(\OWNt;SNRe 

where  is  the  relative  amount  of  focus  (ie.  RJPf)- 

In  order  to  use  (12)  as  the  basis  for  evaluating  the  impact  of  potential  heater  array  errors  on  patch 
location,  each  of  the  error  sources  to  be  evaluated  must  be  tied  to  an  equivalent  additive  noise  error. 
Eor  small  scale  errors,  the  three  potentiai  sources  of  error  can  be  related  to  SNRe  as  shown  in 

Table  2,  In  Table  2,  *)tt)  corresponds  to  the  relative  amplitude  distortion  (i.c.  100  Oyis  the  percent 

rms  distortion)  and  y/(f)  is  the  absolute  phase  error. 


Table  2.  Equivalent  Additive  Element  Noise  £  SNRC‘$  for  Three  Casses  of  Errors1 


Additive  Noise 

iHO 

■Hi. 

jfcf 

4  24 

5 

4 

1 

m 

1 

!  live  SNR f's  have  tweri  derived  by  assuming  dial  1)  ihe  amplitude  Distortion  produces  an  dearie  Odd  intensity 
from  each  element  of  £\( I  *■  e„),  where  e*  is  a  /.cro  mean  error  ami  includes  the  element  phasing  ami  2)  the 
phase  error  produces  an  dement  Odd  intensity  f;Vt/K-02/2)  ♦  *  cJt/H  o^/2)).  where  die  first  term  is  the 

average  value  and  the  second  is  the  random  zero-mean  perturbation  about  the  average.  In  each  case,  the  first  ami 
second  terms  arc  summed  over  all  dements  and  the  mean  and  variance  calculated.  Htc  ratio  of  the  squared  mean  to 
tlic  variance  forms  the  ^suiting  SNR(.  dividing  this  by  die  number  of  etunenu  gives  iltc  equivalait  SNRf. 


11 


2.3.2  Sensitivity  Trade-offs 


In  this  section  we  combine  the  results  of  Table  2  with  equation  (12)  to  obtain  a  family  of  trade-off 
curves  that  indicate  the  expected  rms  AIM  location  error  as  a  function  of  the  error  contributions  of 
the  four  classes  of  error  sources.  For  a  given  error  budget  assigned  to  each  contributing  factor,  the 
total  rms  error  can  be  found  by 

<4  =  V  <4, +  (13) 


where  is  the  error  variation  due  to  factor  i.  Figure  8  shows  the  dependence  of  the  breakdown 
location  error  on  the  element  signal-to-noise  ratio  for  three  different  power  density  gradients. 
Expressing  the  noise  levels  for  t)(/)  and  XO  in  dB  relative  to  the  unperturbed  waveform,  the  impact 
of  these  perturbations  on  Az  can  be  assessed  directly  from  Figure  8.  Figure  9  gives  the  sensitivity 
of  A*  io  the  phase  error,  yKt). 

Table  3  provides  an  assessment  of  the  achievable  performance  levels  using  currently  available 
techno'ogy.  These  numbers  have  been  obtained  from  vendor  specification  data  sheets  and  in- 
housc  experience  v/ith  related  hardware.  The  cumulative  rms  error  from  all  factors  in  Table  3 
provide  a  rms  deviation  in  breakdown  altitude  of  no  more  than  .36  meters.  This  corresponds  to  a 
rms  dev:atiun  across  the  AIM  patch  of  less  than  .06  wavelengths  for  a  30  Mhz  radar  frequency. 
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Figure  9.  Sensitivity  of  Breakdown  Location  to  Heater  Element  Phase  Errors 


Table  3.  Impact  of  Achievable  Heater  Element  Performance  Factors  on  AIM  Breakdown 


Error  Factor 

Achievable  I\rrfonmancc 

Equivalent  SNHe 

Az  rms1 

Quantization  Noise 

10  to  l£  bits 

60  dB 

018  m 

D/A  Noise 

~2 1,11 

63  dB 

.013  m 

System  Clock  jitter 

(§>  350  Mhz 

N/A2 

N/A 

Power  Amplifier  Distortion 

<  "40  dB  THD  <f  3STOT 

*3 

o 

T 

A 

<  ,\1  m 

Phase  Deviations 

<  i8nm 

>  35  dB 

.13  m 

2.4  AIM  Variations  Due  to  Finite  Step  Heater  Scanning 

In  the  following  discussion,  wc  address  the  issue  of  AIM  irregularities  due  to  the  heater's  finite 
step  beam  motion  in  the  direction  of  the  AIM  inclination.  Figure  10  illustrates  the  step  scanning  of 
(he  heater  beam  along  a  short  segment  of  the  AIM  contour.  In  this  figure,  the  heater  dwells  for  a 

given  lime,  t  seconds,  and  then  steps  the  focal  point  a  third  of  a  bcamwidth  in  the  cross  beam 
direction  and  a  corresponding  distance  in  the  bore-sight  direction.  This  dwell  and  step  process  is 
repeated  for  the  extent  of  the  AIM.  As  figure  10  shows,  while  the  desired  constant  electron 
density  contour  is  a  line  inclined  at  45°,  the  actual  contour  due  to  the  quantization  of  the  beam  step 
has  a  ripple  that  oscillates  about  tire  desired  contour,  litis  irregularity  in  the  AIM  will  produce  a 

1  Assumes  breakdown  power  »s  3  dB  below  tire  focus,  4  =  1/2. 

*  System  clock  jitter  will  affect  all  elements  uniformly  and  will  therefore  result  m  a  timing  error.  Ian  not  a  |xiwvr 
density  perturbation. 
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corresponding  ripple  in  the  phase  front  of  the  refracted  radar  wave.  The  extent  of  wavefront  ripple 
will  depend  Both  upon  the  magnitude  of  the  contour  ripple,  as  well  as  the  depth  of  the  plasma 
density.  A  deeper  plasma  provides  more  refraction  and  correspondingly  more  smoothing  of  the 
AIM  irregularity,  thus  less  phase  front  distortion. 


The  objective  of  this  section  is  to  quantify  the  magnitude  of  these  contour  fluctuations  and 
determine  what  beam  step  sizes  are  required  for  acceptable  radar  performance.  These  irregularities 
are  studied  both  from  a  theoretical  perspective,  as  well  as  from  a  straight-forward  numerical 
simulation  of  the  heater  step  scanning  and  consequent  ionization  along  the  desired  constant  electron 
density  contour.  Quantitative  results  are  provided  for  the  analytic  derivations  and  compared  to 
electron  density  irregularities  produced  by  the  numerical  simulation.  These  results  indicate  that, 
within  the  nominal  parameter  regime  for  the  current  AIM  conceptual  design, 

1 )  the  magnitude  of  electron  density  irregularities  falls  off  rapidly  for  paint  step  sizes  less 
than  1/2  beam  width  and 

2)  a  1/4  beam  step  size  produces  a  density  fluctuation  with  magnitude  less  than  .1%  trl  the 
average  and  a  corresponding  location  fluctuation  less  than  .1  in. 

The  next  section  presents  a  derivation  of  the  analytic  relationships.  These  relationships  arc  then 
used  in  Section  2.4.1  to  derive  trade  off  expressions  for  a  specific  heater  beam  example.  Section 
2.4.2  uses  this  result  to  present  a  quantitative  trade-off  between  beam  step  size  and  resulting 
density  fluctuations.  In  addition,  Section  2.4.2  compares  the  theoretical  results  with  a  numeric 
simulation  of  a  stepped  heater  beam  and  associated  ionization. 
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2.4.1  Analysis  Of  Stepped  Beam  Induced  Density  Fluctuations 

In  the  following  discussion,  the  heater  array  is  assumed  to  be  centered  at  (x,z)  =  (0,0)  and  is 
painting  a  patch  by  stepping  the  beam  focus  from  (-  DJI,  zj-  yDJ 2)  to  (DJ2,  yDJ2).  As 

Figure  1 1  illustrates,  each  step  moves  the  beam  uniformly  in  x  and  z,  (&,  &)  =  (A,  yA). 

The  power  density  of  the  heater1  can  be  described  at  the  nth  beam  step  by 

P^x,z)  =  Pfg(x-nA,z-nA^I,  ( 1 4) 

where  g(x,z)  is  the  heater  beam  relative  gain  pattern  (g(0,  z/)  =  1)  in  the  (x,  z)  coordinates  and  Pf  is 
the  heater  power  density  at  the  focus.  For  a  step  size  A  =  jf,  the  cumulative  energy  per  square 
meter  is  found  by  summing  the  contribution  of  each  beam  dwell  at  ( x ,  z) 

Nn 

&x,  z)  =  £  P0g[x  -  nA,  z  -  nAy\  t .  (15) 

a  -  -NI2 

To  simplify  the  following  analysis,  we  assume  that  g(x,  z)  in  the  region  of  the  focus  is  a  function 
of  a  weighted  distance  between  the  point  (x,  z )  and  the  focal  point  (0,  z/), 

gU,  z)  =  g[ax2  +  b{z  -  z/)2) ,  (16) 

where  a  and  b  determine  the  focal  spread  along  the  x  and  z  axes  respectively  (ie.  the  beamwidth 
and  width  of  focus  respectively). 


While  this  form  is  not  completely  general,  it  does  allow  for  a  wide  class  of  beam  patterns,  such  as 

£  _ 

expi-efi)  and  111  the  focal  region  of  the  beam,  where  d  =  'fax2  +  b(z  -  if)2  .  For  a 

pattern  satisfying  (16),  wc  can  show  that 


■  nA,  z 


(17) 


and.  for  a  contour  in  the  (x.  z)  plane  defined  by  (j  -  if)  =  £  +  yx,  wc  have 


l  aJC2  I 

£(*;0  =  4  — 7  +  (o  +  fcy2)  I  nA  - 

\a  +  by  \ 


UX) 


1  In  this  analysis  we  have  ignored  the  effects  of  self  absorption  on  the  heater  wave.  While  this  is  an  tni>>rtam 
phenomenon  during  the  the  patch  creation  process  (leading  to  the  self-limiting  nature  of  the  peak  ciectfci  density ). 
the  effect  on  the  fluctuation  appears  to  be  one  of  limiting  the  electron  density  to  some  maximum,  llur  wni  lead  in 
fluctuations  in  practice  that  arc  less  than  those  predicted  here. 
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If  wc  now  define 


Figure  1 1.  Geometry  of  Beam  Scanning  Effects  Study 


A(u;0  =  ^'"7^2 +  (a  +  *?)  “2|  * 


the  energy  flux  per  square  meter  along  the  contour  can  be  approximated  by 


•0 

-  X  P0linA-x-  -~Y  C  C 


-&oj  '  *  '  ~77y2  ^  X  ^u-nA)du 

-  ti’o  £  c)  exlp-Ka-  x  +  j2x&  ---- 

\  4  Aa  +  brl 


t/v  Mo.c)+  y  2m«/4c)  COJT 1 2^*^  X  ^ 


In  Equation  (20)  H(f;  £)  is  the  Fourier  transform  of  h(u;  £)  in  the  u  variable. 


=  j  Ku;  0  exp(-j2jrfu)  du  .  (21) 


Equation  (20)  shows  that  the  energy  flux  along  contours  parallel  to  the  AIM  surface  consists  of  an 
average  flux,  tPoH(0;  £),  and  a  fluctuation  with  period  equal  to  the  beam  step  size.  The  Fourier 
coefficients  of  the  fluctuation  are  given  by  the  Fourier  transform  of  the  beam  pattern  at  frequencies 

j .  For  beam  patterns  and  step  sizes  of  interest  to  us,  the  majority  of  fluctuation  is  contained  in  the 

first  harmonic,/*  =  1/A  From  (20),  the  first  harmonic  ripple  relative  to  the  average  is 


2H{UA\  0 

m-,o 


(22) 


While  the  energy  flux  represented  by  Equation  (20)  is  related  to  the  ionization  rate  and  consequent 
electron  density,  the  relationship  is  highly  nonlinear  and  therefore  (20)  only  gives  a  sense  of  how 
the  resulting  density  irregularities  will  depend  upon  the  heater  step  size.  Derivation  of  a  more 
accurate  relationship  follows. 

The  electron  ionization  rate  has  been  shown  to  have  the  following  dependence  on  heater  power 
density 


V->/  =  9  x  10 


(1  -G(PIPc)) 


(23) 


G  (PIPC)  =  3.7  x  10-4  exp  [7.9  fFJV). 

where  N  is  the  ambient  neutral  density,  Pc  is  the  critical  power  density  requited  for  breakdown, 
where  the  ionization  equals  rcatiachment  (ie.  (1  •  G(!))  =  0).  Along  the  £  contour  the  power 
density  is  given  by  P0  h(u.  £),  therefore  the  ionization  rate  at  location  x,  £  during  the  nth  twain 

p 

dwell  can  be  written  as  a  function  of  jr  Hu\  £), 

.24) 

Since  the  buildup  of  electron  density  in  the  regime  of  interest  for  AIM  is  of  the  cascade  tyjic,  it  is 

described  by  an  exponential  function  of  Jv,(f)dl.  Using  (24)  for  the  ionization  rate  during  the  nth 

beam  dwell  and  letting  Na(x,  £)  denote  the  ambient  electron  density,  tlic  electron  density  after  the 
complete  heater  scan  is 
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(25) 


N«(x;C)  =A l0(x;QexJ  ^  vfyflnA-x- 

\n  =  -N!2  yc  \  a  +  kf 


tJ 


1) 


where 


-  JVfl(*  0  vM  +  X  2  v(«M:  C)  ^2^  +  M0)|  j 

-  ^(x;  o  «p(t  vfe  d)|  ■ I  +  f  X  2  W/i;  0  cos  j2^x  +  ^Oj 

WiO=j  v^4u;Cl}ejcp(-j2zfu)du. 


(26) 


The  steps  for  deriving  (25)  are  analogous  to  Equation  (20).  Equation  (25)  shows  that  the  electron 
density  along  contours  parallel  to  the  AIM  surface  will  exhibit  fluctuations  due  to  heater  beam  step 

scanning.  These  fluctuations  have  a  period  equal  to  the  beam  step  size,  4,  with  a  relative  first 
harmonic  fluctuation  magnitude  proportional  to  the  beam  dwell  time  and  the  value  of  the  Fourier 
transform  in  (26)  at/pl/A  If  we  let  the  average  electron  density  on  contour  £  be  defined  as 

Ngjivt  ~  Na(x\  0  exp( x  v>(0;  0) ,  (27) 

the  relative  fluctuation  in  the  first  harmonic  is 


SNe 

N€.avt 


=  T2V(l/d;d 

"  wix;  C)J  VM 


(28) 


Equation  (28)  corresponds  to  the  ratio  of  the  peak  electron  density  fluctuation  along  the  desired 
constant  density  contour  to  the  average  density.  In  order  to  quantify  (28).  we  must  specify  a  heater 

beam  pattern,  g(ax2  +  biz-z/)^),  an  ionization  rate  vs (P),  and  then  evaluate  the  Fourier  transform  in 
Equation  (26).  While  this  process  docs  not  explicitly  account  for  the  change  in  beam  pattern  due  to 
absorption,  it  does  quantify  the  density  variation  during  pre-clamping  stages  of  ionization. 


2.4.2  Art  Illustrative  Example 

In  the  following  discussion  we  use  the  results  in  Section  2.4.1  with  a  specific  heater  beam  pattern 
to  develop  explicit  relationships  for  Equations  (22)  and  (28).  For  the  purpose  of  illustration  «c 
assume  a  Gaussian  beam  pattern  in  the  focal  region  of  the  heater.  For  a  3  dB  half-width  of  Bt  and 
B;  in  the  x  and  2  directions  respectively,  the  heater  beam  pattern  is  expressed  as 


g(X.  2) 


l2M) 
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Using  (16)  and  (29)  to  identify  the  appropriate  expressions  for  a  and  b,  (29)  can  be  written  in  the 
same  form  as  (16);  hence,  equation  (19)  gives 


M2LK2  | 

mw  +  *el)  | 

\B}+?B2j 

C UH U 

\  B}  B\ 

/»(«;  0  =  exp^ 

The  Fourier  transform  of  (30)  in  the  u  variable  may  be  performed  anaMically,  producing 

I  Mr 


(30) 


m  o= 


nB}B} 


exa 


\B}  +  ?B}  I 


exp 


n2B}Bl 


n 


(31) 


Substituting  (31)  into  (22)  and  simplifying  gives  the  relative  energy  fluctuation  along  the  AIM 
surface  contour 


2/f(lM;  Q 

mo-.Q 


n2  B}B}  /u2\ 
//«(!-)  B}  +  fBl  U)  ] 


(32) 


2  Q 

f  igure  12  shows  a  plot  of  8e  as  a  function  of  the  number  of  heater  scanning  increments  jicr 

bcatmvidth.  Although  this  plot  is  for  specific  parameter  values,  Bx  =  36  m  and  Bz  =  2,700  m,  the 
results  arc  widely  applicable  across  the  potential  operational  regime  of  the  AIM  heater,  since,  for 
Bx  «  Bx, 


(33) 


As  Figure  12  illustrates,  the  fluctuation  in  ionization  energy  along  a  constant  density  centum 
becomes  extremely  small,  for  even  moderate  beam  overlapping;  for  example,  with  only  2 
increments  per  beamwidth  the  fluctuation  in  energy  is  59  d!J  below  the  average  energy.  As  will  bo 
seen  below,  these  low  energy  fluctuations  translate  into  correspondingly  low  electron  density 
fiuc'uations. 


In  ort‘-.r  to  obtain  a  closed  form  expression  for  the  ionization  rate  and  resulting  electron  density,  *e 

I* 

will  approximate  the  ionization  rate  by  a  Kill  order  polynomial  in  Substituting  (30)  into  (33). 

•  e 

defining  ^  such  that  Pc  -  Pa  g(f\  if  *  ^).  and  us  mg  the  polynomial  approximation. 


£ 

=  N*  A  V  ad 


\exfi 


*Mi) 


f 

x 


Ik- 
exp  — 

t 


B}  *y>B}  B}  j 

kl*tiy{B?  +  rB})  \ 


B}  Bi 


(31 1 
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Increments  per  Beam  width 

Figure  12.  Fluctuation  in  Energy  Flux  Versus  Scanning  Increment 


Substituting  equation  (34)  into  (26  and  performing  the  integration  gives  the  following  Fourier 
transform  of  v: 


As  was  indicated  above,  over  the  parameter  regime  relevant  to  AIM,  Bt  «  B}:  therefore,  (3 r>j 
may  be  approximated  by 
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Given  a  K  order  polynomial  approximation  for  the  ionization  rate,  Equation  (37)  gives  the  relative 

electron  density  fluctuation  along  the  desired  constant  density  contour  located  £  meters  from  the 
heater  focus  if.  As  (37)  demonstrates,  the  fluctuation  is  a  function  of  the  heater  focus  dimensions, 

Bx  and  if*,  the  distance  from  the  focus  to  the  altitude  where  ionization  begins,  and  the  heater 
step  size,  A  The  next  section  will  quantify  these  fluctuations  for  a  typical  set  of  AIM  parameters. 

2.4.3  Quantitative  Results 

A  numerical  simulation  has  been  performed  foi  a  specific  set  of  heater  beam  parameters.  The 
simulation  uses  the  one  dimensional  breakdown  model,  as  defined  in  line  1  of  Equation  (25),  for  a 

fixed  value  of  £  over  a  range  of  x  encompassing  the  patch.  Figure  13  shows  the  results  of  this 
simulation  over  a  short  segment  of  the  patch  for  tne  following  heater  beam  parameters  given  in 
Table  4. 

The  fluctuation  level  measured  for  this  simulation  is  compared  to  the  analytic  predictions  derived  in 
the  previous  section.  In  order  to  quantify  Equation  (38).  a  1 2tT|  order  polynomial  approximation  to 


Table  4.  Heater  Beam  Parameters 


Focus  attitude: 

mmmnmm 

Heater  spot  width  at  focus: 

■gst&iiEKBi 

Heater  spot  focal  width  (in  altitude): 

5.500  to 

AIM  Contour  Distance  from  focus: 

r 

s 

-4.675  m 

i 

-5,500  m 

Heater  step  size: 

A 

variable 

■EM 

55  ptscc* 

v(£  )  was  derived  by  selecting  polynomial  coefficients  to  minimize  the  mean-square  error,  over 

*  c 

the  power  density  range  of  interest.  Figure  14  compares  the  tire  resulting  approximation  with  the 

actu.il  curve  for  Equation  (25)  over  the  range  .1  to  10.  This  includes  tire  range  of ^  values 

•  < 

along  the  £ contour,  for  which  significant  ionization  occurs.  Using  the  resulting  coefficients  and 
substituting  the  relevant  heater  parameters  into  (37)  provides  a  quantitative  expression  for  trailing 
off  the  healer  step  size  against  lire  resulting  level  of  electron  density  fluctuation.  Ibis  relative 
fluctuation  can  be  translated  into  a  corresponding  location  perturbation.  When  the  fluctuations  are 


!  Ttus,  tunc  correspunds  to  an  Kxuratioa  time  required  for  a  300  MHz  heater  «.  uh  an  KRP  of  i  55  dB.  this  urtv  it 
Owner  for  higher  frequency,  higher  ERP  heater  corUtgurationv 
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small  scale,  perturbations  are  approximately  linear.  Therefore,  the  location  perturbation  is  related 
to  the  density  gradient  by 


&  = 


•  8n,  ■  . 


(38) 


For  a  linear  gradient  with  depth  D ,  (38)  is  bounded  by  Sz  <  D  8n Figure  15  shows  the 
quantitative  tradeoff  between  the  heater's  scanning  step  size  and  the  resulting  electron  density 
fluctuation  and  corresponding  altitude  fluctuation  for  the  heater  parameters  given  above.  As  this 
plot  indicates,  the  density  fluctuation  drops  dramatically  as  the  scan  step  size  decreases  below  the 

half  beamwidth  size.  For  a  j  beam  increment,  density  fluctuations  are  less  than  .1%,  with 
corresponding  location  errors  less  than  .1  meter. 


5.000  MO*7 
4.000  MO*7 
3.000  MO*7 
2.000  MO*7 
1.000  MO*7 
"I 

I  010  MO*7 
I  000  MO*7 
9.000  MO** 
9  SQO  MO*** 

9.9T6  MU’* 
9.9?2  MO** 
9  968  MG** 
9.«4  MO** 


Horizontal  Dimension  (m) 


Figure  13.  Numerical  Simulation  Of  Electron  Density  Fluctuations  Due  to  Step  Scanning 
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Figure  14.  Polynomial  Approximation  to  Ionization  Rate 


3.0  AIM  Propagation  Effects 


3.1  Background 

We  begin  our  discussion  of  the  radar  system  performance  with  a  brief  derivation  of  a  radar  range 
equation  appropriate  for  AIM  detection  analysis.  The  AIM  geometry  differs  from  that  of  a 
traditional  radar  by  the  introduction  of  the  reflecting  layer  in  the  propagation  path.  This  layer 
intercepts  the  radar  beam  and  refracts  the  intercepted  wave  toward  the  direction  of  the  target. 
Accounting  for  the  reflecting  area,  losses  during  the  reflection  process,  and  the  two  phases  of 
spreading  loss  (pre-  &  post-reflection)  make  up  the  major  difference  between  an  AIM  range 
equation  and  a  line-of-sight,  monostatic  radar. 

The  relevant  radar  system  and  patch  geometry  parameters  are  defined  as  follows: 


Radar  System: 

Pt  -  average  transmitted  power 
Pr  =  average  received  power 
Gt  =  transmit  antenna  gain 
Ar  =  receive  aperture 

/  =  RF  frequency  (K  =  wavelength) 

B  =  signal  bandwidth 

&s  =  angle  of  surveillance  coverage 
Ts  -  scan  time 

Id  =  Ts  01® j  (dwell  time) 

Lx  =  radar  system  losses 
r  -  total  effective  noise  temperature 
(including  atmospheric  noise) 

a,  =  target  cross  section 

AIM  Gcometry/Losses: 

h  =  AIM  altitude 
R  =  range  from  die  AIM  to  target 

0  =  incident  angle  at  the  AIM 

A,  =  azimuthal  dispersion  off  the  AIM 

/V  ®  elevation  dispersion  off  die  AIM 
tv  =  AIM  width 
I  =  AIM  length 
ip  »  AIM  one-way  absorption 

yp  =  one-way  Faraday  rotation  /  Pm, 
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The  appropriate  expression  for  received  power  can  be  best  understood  by  applying  the  Fresnel- 
Kirchhoff  diffraction  integral  across  the  transmitting  aperture  to  calculate  the  field  at  the  patch  and 
then  across  the  patch  to  calculate  field  reflected  through  the  patch  to  a  target  at  a  distance  R  from  the 

patch  in  direction  of  the  reflection  angle  0(see  Figure  16)[10,  11].  Assuming  a  constant 
illumination  across  the  antenna,  the  resulting  integration  can  be  expressed  as 


F(P)  =  j  (  I  ds^dsp. 


X R\Rz 


Sp  SK 


where 


R\  =  h  +  x2tand-  *lm2e+&-x'?+iyryif., 
2h  2  h  2/i 


(39) 


(40) 


R2  »‘R  -  x2tan$  -  ^tan26  +  +  ^  (41) 

dsx  =  dxjdyj,  and  dsp  -  secddxidyi-  The  approximate  path  lengths  in  (40)  and  (41)  are  used  to 
approximate  the  relative  phase  delays  (Fresnel  approximation),  while  retention  of  the  first  term  is 
adequate  for  accounting  for  the  MR  spreading  loss.  Substitution  of  (40)  and  (41)  into  (39)  gives 


Squaring  (42)  and  multiplying  by  the  transmit  power  density  at  the  antenna  gives  the  power  density 
at  the  target.  Using  the  principle  of  reciprocity,  the  reflected  power  from  the  target  can  be 
calculated  analogously,  where  the  receive  aperture  is  used  in  place  of  the  transmit  aperture. 

The  quadratic  terms  in  the  phase  expression  of  Equation  (42)  divide  the  (w,  l,  6.  h,  R)  parameter 
space  into  numerous  regions,  depending  upon  the  relevance  of  the  different  quadratic  phases  to  the 

integral.  In  general  if  the  quadratic  term,  xlhR/2X(h  +  /?)<:!,  then  that  term  can  be  omitted  front 

the  phase  argument.  Consideration  of  the  relevant  range  of  values  for  (w,  /.  0,  h.  R),  leads  to  the 
.-moderation  of  4  parameter  regimes: 
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Cl:  (l/2?cos*e(h  +  R)  <  A  Mlfth  +  R)  c  l  C2.  (Ikfcos^h  +  R):  ^  (w/2 f(h  +  R)  g  ^ 
2AA/?  *  *  2AM  *  2AM  ‘  ’  2AM 

C3.  (l/2?cos26(h  +  R)<  u (w/2?{h  +  R)  ^  \  Ci.  Wlfcos^h  +  R)  >  {  , (w/2 ?{h  +  R)>  ] 
2 AM  '  '  2AM  ’  2AM  '  ’  2AM 


These  cases  correspond  to  the  radar  antenna  being  in  the  Fresnel  region  (case  4),  far-field  (case  1), 
or  some  combination  thereof  (cases  2  and  3),  relative  to  the  AIM.  For  example,  a  50  MHz  radar, 
with  an  AIM  created  at  70km  and  a  target  at  1000  km,  will  come  into  the  Fresnel  region  when  / 
exceeds  790  m  and  w  exceeds  560  m.  These  cases  are  important,  because  they  determine  the 
appropriate  form  for  the  field  intensity  range  equation,  as  follows: 


Cl:  |F(P)| a- 


%h  +  R) 


C2:  |F(P)|  a 


AJcosd 


X^2ihR(h  +  R) 


C3:  \F(P)\a 


Axw 


AhR(h  +  R) 


C4:  \F(P)\  a 


Axwlcos0 
A  2hR 


In  order  to  achieve  the  azimuth  and  elevation  beam  spreading  consistent  with  the  desired  coverage 
(1.5°  to  3°  and  15°  respectively),  the  AIM  surface  profile  will  be  slightly  curved.  This  curvature 
will  cause  different  portions  of  the  patch  to  reflect  along  different  azimuth-elevation  directions, 
filling  in  the  entire  azimuth-elevation  beam.  For  the  purposes  of  system  performance  analysis,  the 
propagation  can  be  modelled  by  assuming  that  the  AIM  and  radar  antenna  are  both  in  the  far-field 
with  respect  to  each  other  (case  4).  The  resulting  model  can  be  partitioned  into  the  four  segments 
of  the  radar  wave's  propagation  to  and  back  from  the  target.  These  four  segments  arc  transmitter- 
lo-paich,  patch-to-target,  target-to- patch,  and  patch-io-rcceivcr.  Each  corresponding  component  of 
the  range  equation  accounts  for  the  spreading  loss,  absorption,  beam  dispersion,  and  relevant  cross 
sections  along  the  related  propagation  path.  The  power  density  at  the  end  of  each  segment  is 
summarized  by 


f  Pi  Ci 

Ap  Gp 

<7; 

Ap  Gp 

U;r  /i2. 

4  kR1 

AkR2. 

An  k* 

transmit 

-to* 

patch 

pud) 

-u> 

turgci 

target 

•10- 

jwich 

patch 

•U> 

receiver 

(43) 


In  (43)  the  effective  AIM  reflecting  aperture  and  gain,  Ap  and  Gp  respectively,  arc  given  by 


Ap-Lpwl  cos(0) 


(44) 


and 


Gp  =  . 

4^  Af 


(45) 
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For  a  coherent  dwell  and  receiver  noise  given  by  kT,  the  signal-to-noise  ratio  at  the  radar  is 


SNR  = 


PrTd 
kT  ' 


Substituting  equation  (43)  into  (46)  and  combining  terms  gives 


SNR  = 


P,  G,  Aj  Gp  (7t  Ar  Ls  Td 
(4 af  h4  R 4  kr 


(46) 


(47) 


Equation  (43)  assumes  that  any  losses  due  to  Faraday  rotation  during  the  wave's  transit  is 

negligible  (i.e.  yp  «  1).  However,  in  the  case  of  clutter  backscatter  (particularly  sea  clutter),  even 
small  amounts  of  rotation  can  result  in  significant  increases  in  clutter  returns.  This  is  due  to  the 

large  differences  between  horizontal  and  vertical  clutter  cross  sections  (cr/i/,  &  Gw  )•  Typical  ratios 

of  —  are  on  the  order  of  20  dB  or  more.  Taking  these  cross  sections  into  account,  the  signal-to- 
®hh 

clutter  ratio  (prior  to  clutter  cancellation)  for  a  horizontal  transmit  and  horizontal  receive 
configuration  in  sea  clutter  is  given  by 


SCR  = 


cos  (y) 


(48) 


where  B  is  the  radar  waveform  bandwidth,  c  is  the  speed  of  light,  and  yris  the  grazing  angle  at  the 
surface.  As  will  be  discussed  below,  the  small,  predictable  amount  of  polarization  rotation 
imposed  by  the  plasma  layer  allows  an  AIM  system  to  exploit  the  large  difference  between 
horizontal  and  vertical  backscatter. 


3.2  Absorption  and  Dispersion  Losses 
3.2.1  Absorption 

The  absorption  of  the  radar  wave  during  reflection  has  a  major  influence  on  the  design  and 
performance  of  the  radar  system,  and  the  need  to  create  patches  with  low  absorption  strongly 
affects  the  design  of  the  heater  system.  The  absorption  of  a  patch  is  determined  by  the  electron- 

neutral  collision  frequency  vm,  and  the  distance  which  a  ray  travels  within  the  plasma.  This 
distance,  in  turn,  is  determined  by  the  radar  frequency  and  the  electron  density  profile. 

An  intuitive  appreciation  of  the  importance  of  these  factors  in  determining  absorption  may  be 
gained  by  examining  an  analytic  result  for  tltc  power  absorption  in  decibels  of  a  patch  whose 

plasma  density  varies  linearly  with  distance  into  the  patch,  until  (op  =  ci)  (plasma  frequency  equals 
the  radar  frequency)  at  a  depth  £>l. 


1  A  derivation  of  this  result  can  be  found  in  ( 1 2]. 


(49) 
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For  typical  AIM  geometries,  the  angle  of  incidence  (0)will  be  35-45°.  Although  a  linear  plasma 
density  variation  is  not  typical  of  those  expected  to  be  produced  by  AIM  heaters,  this  expression  is 
useful  for  the  qualitative  insight  it  provides.  Since  8  is  fixed  by  the  geometry,  the  absorption  is 
determined  by  D/A  and  vjtir,  if  their  product  is  small,  the  overall  absorption  will  be  small.  For  a 
given  radar  RF,  the  minimum  value  of  vja  is  determined  by  the  ambient  collision  frequency,  and 
D/A  is  determined  by  the  focusing  capabilities  of  the  heater  system. 

Immediately  after  the  patch  is  created,  the  electron-neutral  collision  frequency  vm  is  very  high.  The 

electrons  in  the  patch  rapidly  (on  the  order  of  1  }is)  transfer  their  energy  to  the  neutral  atmosphere, 
and  the  collision  frequency  returns  to  its  ambient  value.  The  ambient  value,  which  is  the  lowest 
possible  collision  frequency  at  a  given  altitude,  is  shown  in  Figure  17. 


Altitude  (km) 

Figure  17.  Ambient  Electron- Neutral  Collision  Frequencies 

The  ambient  collision  frequencies  shown  in  Figure  17  represent  an  important  restriction  on  the 
altitude  of  the  patch  for  practical  AIM  systems.  For  example,  if  we  require  that  vjco  Ire  less  than 
.  1 .  wc  can  sec  that  a  50  MHz  radar  system  is  restricted  to  patch  altitudes  above  60  kilometers. 

The  tintc  required  for  the  patch  to  cool  after  formation  also  varies  with  altitude.  Figure  1 8  below 
shows,  as  a  function  of  altitude,  the  time  required  for  the  electrons  to  cool  (vm  is  within  a  factor  of 
2  of  its  ambient  value).  Because  the  fractional  ionization  is  very  low.  the  cooling  docs  not  result  in 
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any  significant  change  in  the  neutral  temperature.  Note  that  in  terms  of  the  usable  lifetime  of  a 
patch,  which  is  expected  to  be  of  the  order  of  .1  - 1  s,  the  cc  Ming  times  are  quite  short. 

In  addition  to  the  heating  effects  during  ionization,  the  cooled  plasma  may  be  reheated  by  the  radar 
wave  if  its  power  density  is  high  enough.  Since,  the  time  required  for  reheating  is  of  the  order  of  1 

ps  and  decreases  at  higher  powers,  the  use  of  short  radar  pulses  to  avoid  this  effect  is  limited. 

Figure  19  shows  the  ratio  of  vm  after  heating  to  the  ambient  value  for  a  50  MHz  radar  and  a  patch 
altitude  of  70  kilometers. 

Radar  patch  heating  establishes  an  upper  limit  on  the  power  density  which  the  radar  system  places 
on  the  patch.  As  the  radar  power  density  exceeds  the  threshold  level  for  heating  (about  0.1  W/m2 
in  the  example  above),  the  relationship  between  incident  and  reflected  power  becomes  nonlinear. 
At  some  point,  a  maximum  reflection  is  reached,  where  additional  power  input  causes  additional 
absorption  sufficient  to  cancel  the  power  increase.  Above  this,  additional  incident  power  actually 
causes  the  reflected  power  to  decrease.  This  limiting  effect  is  shown  in  Figure  20. 
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Figure  20.  Self-Limiting  of  Reflected  Radar  Power  by  Patch  Heating 
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3.2.2  Dispersion  Loss 

Random  errors  during  the  patch  formation  process  may  result  in  patches  which  have  small-scale 
departures  from  the  desired  shape.  The  causes  of  these  errors  and  their  expected  magnitude  have 
been  discussed  in  Section  2.  In  this  section,  we  examine  the  impact  of  heater  beam  random 
positional  errors  on  the  radar  system. 

In  order  to  assess  the  impact  of  errors  on  the  ultimate  performance  of  the  AIM  radar  system,  we 
first  derive  the  relationship  between  rms  patch  irregularity  and  the  resulting  reduction  in  effective 
patch  gain.  Figure  20  illustrates  the  relevant  two  dimensional  geometry  and  associated  parameters 
for  the  AIM  cloud. 


Using  the  ray  path  model  indicated  in  Figure  21,  the  direction  of  the  exiting  ray,  flU,  2)  is  related 
to  the  spatial  pltase  gradient.  V#*, :),  of  die  propagating  wave  by 


£(*.  *) 


*Vd0r,  2) 

k 


(SO) 


31 


where  k  =  2it/X,  is  the  wave  number  and  0(x,  z )  is  the  relative  phase  of  the  wavefront  at  location 
(x,  z).  The  unit  vector  in  (SO),  which  is  locally  orthogonal  to  the  surface  of  constant  phase,  is 
given  by  the  direction  of  the  ray  leaving  the  plasma. 


sin(2£  +  80 x,  z)) 
cos(2(  +  80x,  z)) 


(51) 


At  any  given  point  in  time,  the  relative  phase  of  the  exiting  wave  at  the  plasma  surface  can  be  found 
by  integrating  the  phase  gradient  projected  along  the  direction  parallel  to  the  plasma  surface.  In 
terms  of  the  parameters  indicated  in  Figure  21,  the  relative  phase  along  the  AIM  surface  is 

~  4^sin(0,  £cos(£))) 

=  <Po  -  k  I  cos(C  +  sin(£),  u  cos(O))  dw  (32) 

K 

»  0o*^(^-^p)cos(O  +  isin(O|  50(usin(£),  ncos(£))du  . 

K 


wliere 

C  =  the  distance  along  the  AIM  surface  (in  the  plane  of  inclination), 

£p  =  the  distance  from  the  point  of  entry  to  the  point  of  exit,  for  the  unperturbed  plasma, 
and 

80 x,  z)  =  the  direction  error  of  the  exiting  wavefront  at  (x,  z). 

For  a  linear  plasma  profile,  D.  Mosher  { 1 3}  has  shown  that 

£p  =  4D  coi(0  (53) 

where  D  is  the  patch  depth  (distance  from  the  ambient  electron  density  to  a  plasma  frequency  equal 
to  the  radar  frequency).  In  the  case  of  the  linear  electron  density  profile,  Mosher  also  shows  that 
die  direction  of  the  exiting  wave  from  is  related  to  the  positional  perturbations  of  the  AIM  profile 
by 


8&1  sm(0.  (*  cos(O)  = - i  - (c(f)  *  t\(  -  Lp)\ .  (54) 

m  com 

where  £(f)  is  the  position  perturbation  relative  to  the  ideal  AIM  profile,  orthogonal  to  the  AIM 
surface.  For  the  purposes  of  our  analysis,  wc  assume  that  c  position  errors  arc  on  the  same  order 
as  Az,  thus  the  approximation  eff)  =  dz(£). 

The  reduction  in  effective  gain  from  the  AIM  surface  due  to  the  wave  front  perturbations  is  given 
by  the  integration  of  the  electric  field  at  the  plasma  surface,  which  can  be  reduced  to 
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exp(j<p(C)  -jk{L  -  f)cos(O)  dC 


(55) 


where  L  is  the  length  of  the  AIM  surface,  the  phase  term  is  compensated  relative  to  the  far-field 
path  length  in  the  direction  of  the  maximum  gain,  and  integration  is  performed  over  the  extent  of 
the  AIM  surface  associated  with  an  exiting  wave  front.  Substituting  equation  (52)  into  (55)  and 
simplifying  gives 


whore 


<5 <p(t)  =  k  sin( 0 1  8&u  sin(C ).  u  cos(Q)  du  . 
Jt * 


(56) 


(57) 


Sn.ce  MO  is  a  function  of  Az(t)  and  Ar  has  been  modelled  as  a  random  perturbation,  MO  is  best 
characterized  through  its  statistical  moments.  The  first  order  moment  is  found  by  taking  the 
expc\  t  i  value  as  follows: 


fluol  -  r 

(t 
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where 


v)  =  £j|*<u)  -  6ftv)f} 

e  L  jl"J9h  f  h  .  Ikjt)  Rs{(u  -  v)  ♦  r))  dx 
2D2coiHQ  I  I?} 


(59) 


In  (59),  Rf(  r)  is  the  autocorrelation  function  of  S(f),  or  in  our  case  of  Ai [(). 
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Therefore,  given  a  statistical  model  for  the  AIM  positional  errors,  the  impact  of  these  errors  can  be 
related  to  a  corresponding  reduction  in  tne  unperturbed  patch  gain.  This  loss  in  gain  can 
subsequently  be  expressed  in  terms  of  its  impact  on  the  overall  radar  detection  performance.  In 

order  to  quantify  the  loss,  we  assume  an  exponential  autocorrelation  for  e(l), 

RAD  =  £|f#+  t)£(£)| 

=  of  exA -.69^/j2) ,  m 

where  of  is  the  rms  location  error  and  Le  is  the  50%  correlation  distance.  Equations  (60)  and  (59) 
can  be  substituted  into  (58)  and  integrated  numerically  to  determine  the  loss  in  patch  gain  as  a 
function  of  the  magnitude  and  correlation  length  of  patch  roughness. 

In  order  to  assess  the  overall  impact  of  AIM  related  losses  on  the  radar  performance,  the  dispersion 
factor,  (58),  and  absorption  factor,  (49),  are  multiplied  to  obtain  a  combined  one-way  AIM  loss 
factor.  In  considering  the  two  components  of  this  combined  loss,  wc  see  that  as  the  patch  depth, 
D,  is  increased  the  dispersion  loss  decreases  (due  to  increased  averaging),  while  the  absorption 
loss  increases.  This  phenomenon  introduces  a  trade-off  between  the  two  competing  loss 
mechanisms,  resulting  rn  an  optimal  AIM  density  gradient  for  a  given  set  of  heater/radar  design 
parameters.  In  comparing  Figure  9  with  Figure  15  it  is  apparent  that  a  heater  step  size  of  1/4 
bcamwidth  will  provide  sufficient  AIM  smoothness,  so  that  the  limiting  factor  determining  beam 
dispersion  will  be  the  heater  control  errors  shown  in  Figure  15.  Figure  22  illustrates  the  loss  trade¬ 
off  between  absorption  and  beam  dispersion  as  a  function  of  patch  depth.  In  this  figure  the 
electron  density  profile  is  a  assumed  to  be  linear  with  altitude  and  the  patch  depth  is  measured  as 
the  distance  from  lire  ambient  density  to  the  plasma  density  corresponding  to  the  critical  density  l  tv. 

(Of,  =  OJ). 


Figure  22.  Typical  Trade-off  For  AIM  1 -os.se  s  vs  Depth  of  AIM  Cloud 
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3.3  AIM  Lifetime  and  Wind  Effects 


3.3.1  AIM  Lifetime 

The  time  of  utility  of  the  AIM  cloud  is  controlled  by  two  physics  considerations.  The  first  of  these 
is  the  cooling  rate  of  the  electrons  following  formation  of  the  cloud.  The  reason  for  this  is  that  the 
fractional  absorption  of  the  radar  signal  is  high  for  energetic  electron  distributions  such  as  required 
for  the  formation  of  AIM  clouds.  Following  the  ionizing  heater  pulse  the  electron  distribution 
relaxes  in  energy  by  first  exciting  optical  lines  and  subsequently  vibrational  and  rotational  levels  of 
N 2.  Cooling  rate  estimates  are  of  the  order  of  10  5  sec  and  are  consistent  with  laboratory 
experiments  [14),  Further  analysis  of  the  cooling  processes  is  planned  using  the  Fokker-Planck 
code.  The  useful  life  of  the  AIM  cloud  starts  following  the  electron  relaxation,  since  at  this  point 
the  fractional  absorption  of  the  radar  signal  is  minimized. 

The  end  of  the  AIM  cloud  useful  lifetime  (the  second  physics  consideration)  is  controlled  by  cither 
local  density  decay  processes  or  by  non-local  transport  process.  For  the  range  of  altitudes  under 
consideration  (60-80  km)  the  dominant  transport  process  is  related  to  neutral  winds.  Electron 
diffusion  is  extremely  weak  since  it  is  ambipolar.  Therefore,  lifetime  estimation  has  concentrated 
on  the  local  decay  processes  due  to  their  shorter  time  scale.  Previous  work  on  the  subject 
( Lombardini,  ( 17)  and  Gurevich,  [16])  identified  the  three-body  attachment  process  as  the  cause 
of  the  electron  density  decay  of  the  electron  distribution  after  it  has  relaxed  in  energy,  which 
involves  two  O2  molecules  and  one  electron.  The  lime  scale  for  this  process  is 


where  is  the  neutral  density  of  O2  molecules  (#/cm~).  Therefore,  depending  on  altitude,  it  is 
of  the  order  of  .01-1  sec.  During  a  review  of  the  literature,  it  was  noted  that  in  several  laboratory 
experiments,  bosh  in  the  US  and  USSR,  in  which  the  decay  of  plasma  due  to  local  processes  v  ns 
observed  (17-211  the  plasma  lifetime  was  much  longer  than  expected  by  the  value  of  fa  given 
above.  Furthermore  the  lifetime  scaled  with  the  square  of  the  plasma  density  square,  indicating 
recombination  as  a  dominant  effect.  Finally  the  lifetime  had  a  dependence  on  the  energy  expended 
tn  producing  the  discharge.  All  of  the  above  indicate  that  detachment  due  to  the  interac’ion  of 
excited  Nj  and  O’  molecules  with  O'  was  balancing  three  body  attachment.  This  resulted  in 
recombination  becoming  the  dominant  plasma  loss  process.  Fur  this  to  occur  sufficient  energy 
should  be  pumped  and  stored  in  vibrational  N;,  thereby  accounting  for  the  dependence  of  the 
plasma  lifetime  on  the  energy  expended.  In  this  case  the  decay  time  will  be  given  by  either 
recombination  or  decay  of  the  vibrational  states  of  Oj  which  are  substantially  longer  (tens  of 
seconds). 


3.3.2  Wind  Effects 

Because  the  AIM  radar  system  relies  on  Doppler  processing  for  clutter  rejection,  the  frequency 
domain  effects  of  the  reflection  process  are  important.  Three  major  causes  of  frequency-domain 
modifications  have  been  identified:  patch  motion,  wind  shear,  and  turbulence.  It  is  important  to 
note  that  all  three  of  these  effects  are  driven  by  motions  of  the  neutral  atmosphere;  because  of  the 
very  low  fractional  toni ration,  I  O'7  or  less,  and  the  relatively  high  neutral  density,  plasma 
oscillations  arc  expected  to  be  very  heavily  damped  and  have  negligible  effect  on  the  reflected 
signal. 
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As  soon  as  it  is  created,  the  ionized  region  will  move  with  the  neutral  atmosphere  in  which  it  is 
embedded.  This  motion  will  cause  a  doppler  shift  in  the  signal  reflected  from  the  patch,  which 
must  be  compensated  by  the  radar  signal  processor.  The  situation  is  similar  to  that  encountered  in 
an  airborne  MTI  radar,  and  signal  processing  techniques  for  these  situations  are  well  developed.  It 
is  not  anticipated  that  motion  of  the  patch  as  a  whole,  regardless  of  the  wind  speed,  will  degrade 
the  performance  of  the  system. 

If  a  wind  shear  exists  at  the  point  where  the  patch  is  created,  the  patch  will  not  only  move,  but  be 
altered  in  orientation  and  shape.  Because  vertical  shear  is  generally  much  more  severe  than 
horizontal  shear,  it  is  expected  to  be  the  dominant  source  of  shear-induced  doppler  effects.  The 
vertical  extent  of  the  patch  is  typically  300-600  meters;  at  altitudes  of  70  kilometers,  vertical  shears 
of  as  much  as  0.05-0.07  sec-1  are  known  to  exist.  (This  figure  represents  the  1%  probability 
maximum  shear  at  Cape  Canaveral,  from  {22].) 

This  worst-case  situation  represents  a  wind  velocity  change  of  15-40  m/s  over  the  vertical  extent  of 
the  patch,  which  corresponds  to  a  rotation  of  0.7-3.6  degrees  per  second.  This  will  cause  a  change 
in  the  direction  of  the  reflected  beam  of  1. 4-7.2  degrees  per  second.  Although  this  has  minor 
effect  on  the  pointing  of  the  beam,  it  does  impose  a  scanning  modulation  on  the  target  signal  given 
by 


<7/  =  0.265^.  (62) 

1  <P 

where  Of  is  the  frequency  deviation  of  the  power  spectrum  in  Hz,  <pKiSX  is  the  scanning  rate  in 
degrees  per  second,  and  0  is  the  bcamwidth  in  degrees.  The  maximum  frequency  spread  for  the 

example  system  discussed  above,  for  which  <plcaj1  =  12°,  is  0.16  Hz.  If  the  coherent  integration 
time  for  the  radar  system  is  less  than  6  seconds,  the  frequency  resolution  of  the  system  will  be  too 
coarse  to  detect  this  modulation.  For  a  radar  frequency,/,  this  frequency  deviation  corresponds  to 
a  velocity  deviation  of 


ov 


C  Of 

=  ~2f  ' 


(63) 


which  for  our  example  equals  0.48  meters  per  second.  This  is  the  worst-case  frequency  spread 
due  to  shear  for  the  example  system;  because  it  is  equal  to  tire  velocity  spread  of  a  calm  sea.  it 
would  have  a  slight  effect  on  clutter  rejection  for  targets  over  sea.  In  general,  the  impact  of  shear 
induced  scanning  modulation  is  expected  to  be  slight. 


3.4  Faraday  Rotation 

Depending  on  the  geographical  location  and  orientation  of  the  patch,  the  geomagnetic  field  will 
cause  some  amount  of  Faraday  rotation  of  the  polarization  of  the  reflected  radar  signal.  In  some 
cases,  such  as  a  system  viewing  targets  over  land,  this  may  be  of  little  importance.  As  discussed 
in  Section  3.1  scattering  characteristics  of  sea  clutter  are  strongly  polarization  dependent.  Figure 

23  shows  the  HF  vertical  and  horizontal  co-polarization  cross  sections  (otv  and  Okh  respectively) 
for  a  moderate  level  sea  state,  based  upon  a  widely  accepted  model  for  back  scatter  from  the  sea 
(23].  These  curves  show  to  be  on  the  order  ut  20  dB  down  from  Ovv.  Therefore,  the 
preferred  polarization  for  reducing  sea  clutter  returns  is  a  linear  horizontal  electric  field.  From 
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Figure  23.  Sea  Clutter  Cross  Sections  for  a  Moderate  Sea  State1 


Equation  (48)  we  see  that  as  long  as  the  Faraday  rotation  can  be  kept  to  less  than  10%  (ie.  yp<  .  1 ) 
the  horizontal  backscatter  will  be  the  dominant  clutter  return.  This  assessment  assumes  that  the 
polarization  separation  is  limited  by  the  rotation  in  the  AIM,  rather  than  at  the  transmit  and  receive 
channels. 

In  modeling  Faraday  rotation,  it  is  necessary  to  employ  a  vector  wave  equation  to  detennine  the 
change  in  the  polarization  of  the  incident  wave.  The  approach  selected  follows  that  of  Budden, 
[24],  i.e.  the  calculation  of  the  election  coefficient  rather  than  the  electric  field  vector  to  avoid 
singularities  at  the  reflection  point. 

In  this  approach,  a  vacuum  reflection  coefficient  matrix  R  is  defined  to  be 

[&l  ■  * 

Where  Eh  and  Ev  are  the  equivalent  vacuum  fields  at  the  point  where  R  is  evaluated,  and  the  i  and  r 

subscripts  indicate  incident  and  reflected  waves.  The  R  matrix  indicates  the  degree  to  which 
horizontally  polarized  incident  waves  are  transfoimcd  into  vertically  polarized  reflected  waves,  and 
vertically  polarized  incident  arc  transformed  into  horizontally  polarized  reflected,  by  Faraday 

rotation.  A  complex  matrix  differential  equation  for  R  in  terms  of  altitude  may  be  derived  (22), 

2;'-~=  W2\  +  W22R  ~RW\\  -RW\2R  .  (61) 

dz 

where  the  w's  are  matrix  functions  of  the  plasma  properties  at  that  altitude  and  the  intensity  and 
orientation  of  the  magnetic  field.  This  set  of  equations  may  be  solved  using  any  of  the  various 


1  Cross  section  models  taken  Iron)  (23 1 . 
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approaches  for  ordinary  differential  equations;  in  this  case,  an  adaptive  altitude  step  Runge-Kutta 
method  was  used. 

Computations  using  typical  profiles  produced  by  simulations  of  the  plasma  formation  process 
show  rotations  of  as  much  as  2-3  radians  for  the  worst-case  location  and  patch  orientation.  In 
these  cases,  it  may  be  necessary  to  precompensate  for  Faraday  rotation  by  adjusting  the 
polarization  of  the  transmitted  signal.  Because  the  amount  of  rotation  depends  only  on  the  depth, 
location,  and  orientation  of  the  patch,  the  amount  of  compensation  required  can  be  computed  or 
determined  by  experiment  and  applied  each  time  a  patch  is  created  It  is  important  to  note  that  any 
AIM  radar  system  which  scans  in  azimuth,  and  which  uses  polarization  control  to  take  advantage 
of  the  difference  in  sea  reflectivity,  must  already  have  the  ability  to  transmit  varying  polarizations; 
so  using  this  capability  to  compensate  for  Faraday  rotation  does  not  add  additional  complexity  to 
the  system. 
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4,0  Projected  Threat  and  AIM  System  Performance 


4.1  The  Cruise  Missile  Threat 

In  order  to  quantify  the  potential  performance  of  an  AIM  based  system,  a  specific  surveillance 
mission  and  associated  threat  was  selected.  Because  of  the  potential  use  of  this  type  of  system  as 
an  adjunct  to  a  traditional  skywave  OTH  radar  and  the  importance  of  the  mission,  we  have  selected 
the  continental  United  States  strategic  defense  as  the  candidate  mission,  and  the  corresponding 
threat  is  a  low-flying  air  or  sub  launched  cruise  missile.  Selection  of  this  mission  and  threat 
presents  a  number  of  challenging  problems. 

First,  the  surveillance  environment  consists  primarily  of  long-range  (beyond-the-horizon)  air  space 
over  the  sea.  In  order  to  provide  OTH  coverage,  the  AIM  will  be  created  at  approximately  70  km 
altitude,  giving  a  down  looking  beam.  The  corresponding  grazing  angle  will  vary  from  0*  at  the 
AIM  horizon  (~  1 100  km)  to  15*  at  200  km.  Therefore,  as  Figure  23  illustrates  sea  clutter 
backscatter  will  become  a  limiting  factor  as  the  target  closes  in  range.  In  order  to  help  reduce  the 
sea  clutter,  a  horizontal  polarization  is  preferred  (see  discussion  in  Section  3.4). 

The  second  issue  to  be  considered  is  the  relative  geometry  of  a  low-flying  target  over  a  conducting 
surface,  such  as  the  sea.  As  is  well  known,  this  geometry  can  give  rise  to  severe  multipath  nulls. 
An  effective  approach  to  dealing  with  multipath  is  the  use  of  frequency  diversity.  As  the  transmit 
frequency  is  changed,  the  multipath  nulls  move  in  range  (with  the  exception  of  the  first  null  at  the 
AIM  horizon).  This  effect  is  illustrated  in  Figure  24,  where  a  single  frequency  multipath  response 
is  compared  to  the  average  response  due  to  a  transmit  frequency  that  is  hopped  between  three 
frequencies:  43  MHz,  49  MHz,  and  54  MHz.  As  these  curves  illustrate,  the  use  of  frequency 
diversity  all  but  eliminates  the  multipath  nulls. 

Finally,  the  radar  cross  section  characteristics  of  a  cruise  missile  play  an  important  role  in  how  well 
the  radar  will  be  able  to  detect  and  maintain  a  track  on  the  target.  Figure  25  shows  a  numerically 
calculated  radar  cross  section  of  a  cruise  missile  like  target.  The  target  consists  of  a  metallic  cigar- 
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Figure  25.  Numerical  Model  for  a  Cruise  Missile  Like  Target  Illustrates  the  Advantages  of  Having 

Positive  Frequency  and  Polarization  Control1 


shaped  cylinder  with  metal  fins.  The  dimensions  of  the  target  approximate  those  of  a  typical  cruise 
missile.  While  the  resulting  cross  section  calculations  do  not  exactly  match  those  of  a  real  cruise 
missile,  they  do  give  a  good  indication  of  the  relative  dependence  upon  frequency  and  polarization. 
As  these  curves  indicate,  resonance  occurs  in  the  50  -  100  MHz  range,  with  a  consequent  fall  off 
as  f  below  this  region  (corresponding  to  Rayleigh  scattering).  Additionally,  for  low  grazing 
angle,  dosing  geometries  there  is  a  strong  dependence  on  polarization,  with  horizontal  polarization 
providing  from  10  to  20  dB  advantage.  Therefore,  the  positive  frequency  and  polarization  control 
offered  by  AIM  can  be  used  to  simultaneously  exploit  the  cross  section  characteristics  of  both  the 
target  and  the  sea  clutter. 


4.2  A  Baseline  System  and  its  Predicted  Performance 


Baseline  Heater: 

litis  section  describes  the  current  baseline  licatcr  system  and  describes  trade  studies  that  have  been 
performed  to  arrive  at  this  baseline.  There  arc  two  rather  distinctive  approaches  for  simultaneously 
obtaining  the  high  ERP,  achieving  the  positive  field  gradient  needed  for  controlled  breakdown,  anil 
controlling  the  plasma  cloud  with  sufficient  precision  that  it  may  be  used  for  tire  radar  application: 

1 )  a  large,  solid-state  phased  array,  or 

2)  high-power  tubes  (klystrons,  likely)  feeding  several  dishes. 


1  Taken  from  C.  C.  Ou.  Syracuse  Research  Corporation  (251. 
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Our  current  understanding  of  the  phenomenology  and  the  available  technology  leads  us  to  favor  the 
solid-state  approach  for  the  Full  Scale  Engineering  Development  system  because  of  its  greater 
controllability,  likely  lower  life  cycle  costs,  and  greater  reliability.  The  tube  approach  is  currently 
favored  for  an  atmospheric  proof  of  principle  demonstration  because  of  its  lower  development 
cost. 

Based  on  our  study  of  the  plasma  creation  process  and  associated  control  requirements,  formation 
of  an  AIM  that  will  efficiently  reflect  HF/VHF  radar  waves  requires  a  heater  with  baseline  design 
characteristics  as  follows: 

•  Heater  frequency  in  the  UHF  band, 

•  Moderate  (about  10%)  frequency  agility, 

»  Effective  radiated  power  of  1 56  dB W, 

•  Average  radiated  power  of  2-5  MW, 

•  Duty  cycle  of  0. 1-5  %,  and 

•  Flexible  and  accurate  control  of  the  E  field  focus. 

The  current  baseline  concept  for  the  heater  antenna  consists  of  a  large,  phased  array  (about  4  km2) 
with  solid  state  transmit  elements,  each  radiating  at  a  nominal  frequency  of  425  MHz  (Figure  26). 
While  the  number  and  placement  of  the  antenna  elements  is  still  subject  to  a  detailed  design,  the 
general  trade-off  is  one  of  increasing  the  number  of  elements  (ie.  radiating  aperture)  in  order  to 
reduce  the  required  heater  energy  per  AIM  cloud.  As  the  array  area  is  filled,  it  becomes  more 
efficient  in  terms  of  placing  radiated  power  at  the  heater  focus,  consequently  requiring  less  total 
radiated  power. 

To  suppon  the  requirement  for  near-field  focussing  of  the  heater,  each  sub-panel  will  be  digitally 
calibrated  and  synchronized.  Reference  signals  will  be  measured  by  receivers  used  to  provide 
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Figure  26.  lire  FSED  AIM  Heater  Phased  Array  Antenna 


calibration  for  each  sub-panel.  Digital  control  and  waveform  synthesis  will  provide  very  accurate 
synchronization  of  the  transmit  waveforms,  allowing  complete  flexibility  in  synthesizing  the  field 
distribution  at  the  breakdown  region  by  specifying  die  phase  and  amplitude  weights  for  each  sub¬ 
panel. 

An  alternate  approach  would  consist  of  an  array  of  reflectors  up  to  25  meters  in  diameter  driven  by 
high  powered  L  or  S  band  klystrons.  The  number  of  reflectors  may  vary  from  25  to  several 
thousand  depending  on  the  size  of  the  reflector,  AIM  requirements  and  output  characteristics  of  the 
klystron.  The  number  of  klystrons  may  vary  from  10  to  several  hundred  depending  on  their 
output.  A  design  for  a  L  band  klystron  with  1  GW  output  power  for  a  1  ps  pulse  has  been 

proposed.  There  exists  a  commercially  available  L  band  klystron  at  250  MW  with  a  1  ps  pulse. 
This  array  of  reflectors  would  be  focused  as  the  baseline  concept  with  digital  waveform  synthesis 
and  control.  The  ultimate  choice  will  be  determined  based  on  system  trade-offs  of  AIM  creation 
phenomenology,  performance  requirements,  and  cost 


Baseline  Radar: 

The  baseline  AIM  radar  consists  of  a  phased  array  aperture  operating  in  the  30  -  60  MHz  range, 
with  an  ERP  of  98  to  104  dBW.  The  radar  will  operate  with  a  frequency-swept  FM  pulsed 
waveform,  where  the  pulse  width  is  chosen  to  match  the  540  km  two-way  propagation  time,  from 
the  radar  to  the  initial  range  of  surveillance.  This  will  provide  maximum  average  power, 
compressed  range  resolution,  and  co-location  of  the  transmit  and  receive  aperture.  Table  5 
summarizes  the  pertinent  radar/patch  parameters  for  25  MHz  and  50  MHz  baseline  designs.  These 
designs  will  provide  a  360°  azimuth  scan,  with  a  1.5°  to  3°  azimuth  beamwidth  and  range  coverage 
of  200  to  1000  km  for  a  low  flying  target  (50  m  altitude).  Scan  update  times  for  360°  coverage  will 
be  on  the  orticr  of  1  to  4  minutes,  depending  upon  AIM  lifetimes  and  azimuth  beam  width. 


Table  5.  Two  Baseline  Radar  Configurations 


System  Parameters 

IIF  AIM 

VHFAIM 

Average  Power 

Mj 

1.2  MW  H 

Antenna  Aperture  (Cain) 

250  x  250  m  (37  dB) 

250  x  250  m  (43~dB) 

25  MHz 

50  MHz 

1 61kHz 

100  kHz 

Effective  Noise  Temperature 

System  Losses 

10  dB 

lOdB 

Patch  Parameters 

Size 

1 .5  x  3.0  km 

1 .4  x  2.0  km 

Altitude 

76  km 

70  km 

Inclination 

~~  40*  to  45s 

40^0450 

4  dB  (one-way) 

5.5  dB  (one-way) 

Relative  Polarization  (W/ltH) 

-24  dB 

-28  dB 

Beam  Dispersion 

l3razn5Tci 

Baseline  Performance: 

Figure  27  indicates  the  expected  performance  of  the  two  baseline  radars  shown  in  Table  5.  using 
AIM  against  targets  at  200  -  1,000  km  ranges.  These  curves  illustrate  the  sensitivity  of  the  radar 
detection  performance  in  a  moderate  level  sea  clutter  and  in  noise.  All  AIM  and  system  related 
losses  have  been  included  in  the  analysis.  The  minimum  detectable  target  is  defined  as  the 
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minimum  radar  cross  section  required  for  a  probability  of  detection  exceeding  90%  with  a  false 
alarm  rate  of  10"3  (~  100  pre-tracking  false  alarms  per  scan).  As  can  be  seen,  the  VHF  radar  has 
over  30  dB  margin  against  the  target  illustrated  in  Figure  25. 


Target  Range  (km) 

Figure  27.  AIM  Detection  Performance  Provides  30  dB  Margin  at  VHF 
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5.0  Summary 


As  the  results  in  this  paper  have  indicated,  the  operational  geometry  and  radar  characteristics  of  an 
AIM  based  surveillance  radar  can  provide  system  performance  that  is  well  suited  to  the  long-range 
detection  and  tracking  of  small,  low-flying  targets,  as  exemplified  by  the  ALCM  and  SLCM 
threats.  The  performance  characteristics  of  a  stand-alone  AIM  system  are  highlighted  by 

•  200- 1200  km  detection/tracking  ranges  (depending  on  target  altitude), 

•  azimuth  coverage  over  360°,  or  any  fractional  sector, 

•  beam  grazing  angles  of  less  than  10  degrees  (for  65-80  km  high  AIM), 

•  good  horizontal  polarization  control, 

•  operation  at  RFs  in  the  HF  to  low  VHF  (to  about  90  MHz),  and 

•  90%  detection  probability  of  a  -25  dBsm  target  at  1 ,000  km  range,  providing  30  dB 
of  margin  for  typical  low  observable  threat  projections  at  VHF. 

In  the  case  of  strategic  defense  of  the  continental  United  States,  an  AIM  system  can  provide  an 
important  augmentation  to  current  and  planned  deployment  of  more  traditional  skywave  over-the- 
horizon  radars.  AIM  enhances  the  overall  performance  of  these  OTH  radars  by 

•  filling  in  the  range  hole  that  exists  out  to  about  1000  km  due  to  the  minimum  HF 
hop  distance, 

•  mitigation  of  auroral  effects  in  polar  directed  surveillance  sectors, 

•  sustained  operation  through  periods  of  increased  sunspot  activity  and  other  iono¬ 
spheric  degradations, 

•  availability  of  the  upper  end  of  the  HF  spectrum  during  the  diurnal  ionospheric 
cycle,  and 

•  improved  detection  of  LX)  targets  through  frequency  selection  and  positive  polar¬ 
ization  control. 

The  basic  questions  concerning  the  feasibility  of  the  AIM  concept  all  relate  to  how  reliably  one  can 
create  the  AIM  within  tolerances  necessary  for  useful  reflection  of  the  radar  signal.  Issues  that 
directly  impact  AIM  utility  (size,  shape,  orientation,  uniformity,  smoothness,  peak  electron 
density,  steepness  of  the  density  gradient,  and  density  lifetime)  have  been  addressed  above. 
While  some  uncertainties  still  remain,  efforts  to-date  have  reduced  these  uncertainties  to  the  extent 
that  there  appears  to  be  no  phenomenological  issue  preventing  realization  of  the  AIM  concept.  This 
is  bused  upon  extensive  physics  studies  and  system  trade-off  considerations.  Results  of  these 
investigations  indicate  that  the  AIM  can  be  created  using  current  technology,  with  the  necessary 
radar  reflective  characteristics  to  provide  wide  area  surveillance  and  early  detection  and  tracking  of 
the  ALCM/SLCM  threat  of  the  future.  Based  on  these  results,  the  next  phase  of  effort  should 
include  the  careful  design  and  planning  of  an  atmospheric  demonstration  of  the  AIM  concept 
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